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Abstract. We consider a family of open sets which shrinks with respect to 
an appropriate parameter e to a graph. Under the additional assumption that 
the vertex neighbourhoods are small we show that the appropriately shifted 
Dirichlet spectrum of Mg converges to the spectrum of the (differential) Lapla- 
cian on the graph with Dirichlet boundary conditions at the vertices, i.e., a 
graph operator without coupling between different edges. The smallness is ex- 
pressed by a lower bound on the first eigenvalue of a mixed eigenvalue problem 
on the vertex neighbourhood. The lower bound is given by the first transversal 
mode of the edge neighbourhood. We also allow curved edges and show that 
all bounded eigenvalues converge to the spectrum of a Laplacian acting on the 
edge with an additional potential coming from the curvature. 



1. Introduction 

Graph models of quantum systems can often be used to describe in a simple 
way some important aspects of the behaviour of a quantum system. Although 
such models are simple enough to be solvable (because they are essentially 1- 
dimensional) they still have enough structure to model real systems. Ruedenberg 
and Scherr |RuS53j used this idea to calculate spectra of aromatic carbohydrate 
molecules. Nowadays the rapid technical progress allows to fabricate structures 
of electronic devices where quantum effects play a dominant role. Graph models 
like quantum graphs (also called metric graphs) can often be viewed as a good 
approximation of such structures. From the m athematical point of view these 
models were analysed first thoroughly in [ES89J , for recent developments, bibli- 
ography and further applications see [DEOS^ , [KoS99j , |Ku02j or |Ku04j ; note that 
|KaP88j also calculated the eigenvalue asymptotic of a tubular ^-neighbourhood 
of a curve. 

A quantum (or metric) graph is a graph where we associate a length to each 
edge. A natural operator acting on such graphs is given by a self-adjoint extension 
of —(P/dx^ on each edge. We will call such a self-adjoint extension a Laplacian on 
the (quantum) graph. Note that the Laplacian on a discrete graph is a difference 
operator on i2{K) rather than a differential operator acting in 0^. L2{ej). Here, 
K labels the vertices and J the edges Cj, j e J, of the graph. A detailed overview 
on this wide field can be found in |Ku04j or |Ko^ 
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A natural question is in what mathematical sense a quantum graph Mq can be 
approximated by a more smooth space M^. One is interested what Laplacians 
on Mq occur as limit operators from operators on M^. More significantly, 
could be the ^-neighbourhood of an embedded graph Mq C M" or a manifold 
shrinking to Mq as e ^ 0. We call such approximating spaces branched quantum 
wave guides. Recently, spectral convergence in the case of a bounded open set 
Me with Neumann boundary condition has been established in |RSc01j , |KuZ01j 
and |KuZ03j : for an approximation by manifolds see |EP04j . All these examples 
have in common, that the lowest eigenmode of the transversal direction is with 
constant eigenfunction. In this case, the limit operator is the Laplacian on the 
graph with Kirchhoff boundary conditions, i.e., a function / in the domain of the 
Kirchhoff Laplacian is continuous at each vertex and satisfies 

= 0, keK. (1.1) 

In addition, the spectral convergence holds independently of a given embedding 
of the graph. In particular, the convergence is independent of the curvature of 
the embedded edges. 

The case of an approximation by Dirichlet Laplacians on an open set was 
first treated heuristically in |RuS53j . This case is harder to analyse since the first 
transversal eigenvalue equals (F^) = \i/e^ (Ai > 0), i.e., it is of the order 
if e denotes the radius of the cross section = {—s,e) of the approximating set 
M^. A rescaling is necessary, and first order terms of the metric (cf. (j4.3|) ) like 
the curvature become important. In particular, the curvature of the (embedded) 
edge enters in the limit operator as an additional potential. 

Main result. Assume that Mq C is a finite graph. Our aim in this note is to 
show the spectral convergence of the Dirichlet Laplacian on an approximating 
open set D Mq. We suppose that can be decomposed into neighbour- 
hoods Uej of the edges Cj and into neighbourhoods Ve,k of the vertices Vk of Mq 
(cf. Figure m). We assume that Ve^k is e-homothetic to a fixed set Vk. The precise 
definition will be given in Section |21 and El Our basic assumption is that the 
vertex neighbourhoods Vs^k are small, i.e., that 

X^^'iVk) > A?(F) = Ai (1.2) 

where A™(Vfc) is the lowest eigenvalue of the Laplacian A™ of Vk with Dirichlet 
boundary conditions on c^oVfc (i.e., on the boundary induced from the original 
boundary of Mg) and Neumann boundary conditions on djVk, j G Jk, (i.e., on 
the parts where the adjacent edge neighbourhoods labeled by j G Jk emanate, 
cf. Figure I2I). Furthermore, F = (—1,1) and therefore Ai = 7r^/4. We comment 
on this condition in Sectional 
Our main result is 
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Figure 1 . Decomposition of the graph neighbourhood (grey) 
of the graph Mq into edge and vertex neighbourhoods Uej and V^^k- 




V N 

Figure 2. The scaled vertex neighbourhood Vk with the boundary 
part doVk coming from the original boundary and the boundary 
part djVk where the edge Cj emanates. 

Theorem 1.1. Suppose that is an open neighbourhood of a finite graph 
Mo C satisfying the smallness assumption (II. 2j) on each vertex neighbour- 
hood. Denote by Xk{£) the k-th eigenvalue of the Dirichlet-Laplacian /S^^ > 
(counted with respect to multiplicity). Then 

Afc(£) - ^ Afe(O), 5^0 (1.3) 
where Afc(O) denotes the k-th eigenvalue of 
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with A^, = —d'^/dx'j being the Dirichlet Laplacian on the edge Cj = Ij = {0,ij) 
and Kj being the curvature of the embedded edge Cj C (cf. fl4.2|) ). 

Note that the smallness assumption at the junctions Vs^k iniphes that the hmit 
operator decouples, i.e., the hmit operator is the direct sum of operators acting 
on a single edge. In the case of the Neumann Laplacian on decoupling occurs 
if the area of the edge neighbourhood decays faster than the area of the vertex 
neighbourhood; e.g., if the latter scales in each direction of the order with 
< a < 1/2; the vertex neighbourhoods are large obstacles seen from the edge 
neighbourhoods (cf. |KuZ03j or EP04J). In the case of Dirichlet boundary con- 
ditions, in contrast, decoupling already occurs when the vertex neighbourhoods 
scale with e, i.e., even when the edge neighbourhood volume (which is of order 
e) decays slower than the vertex neighbourhood volume (of order e^). 

We also show in Section El that the usual e-neighbourhood := {z G 
I d{z, Mo) < b} does not satisfy our hypothesis since the leading order of 
the lowest eigenvalue is at most fi/e"^ with /i < Ai. Therefore, A^^^ — Ai/e^ has 
a negative eigenvalue tending to — oo (cf. also Lemma 12. 1|) and the conclusion 
of Theorem 11.11 fails. In particular, there is no limit operator on the graph (us- 
ing the simple shift A — * A — Xi/e'^), and the suggestion in jRuS53j . that the 
limit operator is the Kirchhoff Laplacian on Mq, is false (cf. also |Ku02| Sec. 2.1 
and 3.2]). Note that Ruedenberg and Scherr implicitly assumed that the lowest 
eigenfunction does not concentrate around the vertex which is the case as we will 
see in the last section. 

The spectral convergence of a single curved quantum wave guide has already 
been shown in |DE95j and |KaP88j using perturbation methods. Our proof only 



uses variational methods and is a simple adaption of (EP04| . |Ku/OH fKuZ03j or 
|RSc01j . where one compares Rayleigh quotients. 

The paper is structured as follows: In the next section we define the required 
spaces and operators in the case of straight edges. Section |21 is devoted to the 
proof of Theorem 1 1.1 1 in this case. In Section|3]we provide the necessary changes in 
order to prove the result with curved edges. Section El contains some explantation 
on the smallness condition ()1.2|) . and examples where this condition holds or fails. 



2. Preliminaries 

In this section we define the limit space and the approximating space together 
with the associated operators. We first consider straight edges without curvature, 
= 0. In Section m we also allow curved edges. 

Definition of the limit space. Let Mq be a finite connected graph with (metric) 
edges Cj, j E J {cj isometric to an open interval Ij = {0,ij)) and vertices Vk, 
k & K. We denote the set of all j & J such that Cj emanates from the vertex Vk 
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by Jk- The Hilbert space associated to such a graph is 

7i:=L2(Mo) = 0L2(/,) 

which consists of all functions / with finite norm 

ii/iio = ii/iiLo = Eii/^-iil = E / \fA^)\"d^- 

Definition of the limit operator. We define the limit operator Qq via the quadratic 
form 

?o(/):=Ell/ill|=E / \fM)?d^ 
j€J jeJ ''^0 

for functions / G C^(Mo) = @jC^{Ij) (with compact support). The form 
closure of go (also denoted by go) is the extension of go to the closure of the space 
of all such functions in the norm 

ll/llo,i:= ll/llo + ?o(/) 

(see jK66| Chapter VI], |R,S8nj or |Da96j for details on quadratic forms). Note 
that 

domgo = 07ii(/,). 

Remember that fC'il) is the closure of C^°°(/) w.r.t. the norm (||/f + H/'f)^/^. 
The associated self-adjoint, non-negative operator Qo is given by 

Qo = 0A°. (2.1) 

where > denotes the self-adjoint operator —(P/dx^ on Ij with Dirichlet 
boundary conditions. The spectrum of Qq is purely discrete and will be denoted 
by Afc(O), written in ascending order and repeated according to multiplicity. 

Definition of the approximating space. We now describe the family of open sets 
(Me)^, < e < £0; approximating the graph Mo as e — > 0. For convenience 
only, suppose that Mq is embedded in (an abstract definition of in the 
general case will be given soon). Assume that we can decompose into open 
sets containing those points x G Cj with d{x,dej) > aje/2 for some real 
number aj < I/eq and V^r^k 3 ^/c such that the union of their closures equals 
Mg. Here, the edge neighbourhood Usj is isometric to I^j x (both equipped 
with the Euclidean metric) where I^j '■= (0, (1 — aje)£j) and = {—e,e) is the 
scaled cross section. Furthermore, we assume that the vertex neighbourhood ^ 
is e-homothetic to a fixed open set V^. Using a simple coordinate transform we 
have therefore the isometrics 

{U,j,g,^,,) = {ljx F,g,) and {V,,k, 9end) = {Vk, gs) (2.2) 



6 



OLAF POST 



where 



= {I — eaj^dx^ + e^dy^ and g^. 




(2.3) 



are the metrics on the edge resp. vertex neighbourhood. Here, F = (—1, 1) and 
g is the Euchdean metric on V^. In the sequel we use this change of coordinate 
transform without mentioning. Note that the slightly shortened edge neighbour- 
hood is necessary in order to have an embedding for the edge and the vertex 
neighbourhood. 

Although we are mainly interested in the embedded situation as described 
above, we prefer the following abstract setting in order to keep the notation 
of , EP04 and recognise the important geometric objects (not depending on any 
embedding). For each j G J we let Uej be the Riemannian manifold {Ij x F,g^) 
where is given as in ()2.3|) . Here F is the interior of a compact, connected Tri- 
dimensional manifold (m > 1) with metric denoted by dy"^ having purely discrete 
Dirichlet spectrum with first eigenvalue Ai > 0. 

Furthermore, we denote by V^^k the Riemannian manifold {Vk, ge) with g^ = e^g 
where g is a. metric on V^. We assume that 



i.e., the boundary of Vk has as many boundary parts djVk isometric to F as 



Furthermore, we assume that the metric on Vk has product structure g = dx^+dy"^ 
near djVk- 

We can define an abstract manifold by identifying the appropriate boundary 
parts according to the graph Mq. Note that a smooth structure on Mg and also 
a smooth metric g^ of the form ()2.3|) in the respective charts exist since is 
diffeomorphic to a product (0, 1) x F in a neighbourhood of each djVk on both 
sides of djVk, i.e., on Vk and Uj. Strictly speaking we should introduce another 
chart for each j G Jk and k & K covering djVk in order to define the smooth 
structure properly. But since we only use integrals over M^, a cover up to sets 
of measure is enough. The resulting manifold has dimension d = m -\- 1. Note 
that Me need not to be embedded in any space, but the embedded case described 
above is also covered by this setting. 

The associated Hilbert space is 




(2.4) 




L2(M,) = L2(f/.,,) © L2{V,k) 
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which consists of all functions u with finite norm 

|2 



HI? = IHIk = ElHlL + E 

j€J keK 



= ^ / \u\'^{l - aje)e'^ dxdy + / \u\'^e'^dz 

where dy and dz represent the natural measures on F and Vk, respectively. 

Definition of the operator on the manifold. The operator on the thickened space 
we are considering will be the Dirichlet Laplacian on M^, i.e., = > 0. 
The corresponding quadratic form /i^ is given by 

= [ ^—^\d^u\^ + \\dyu\^]{l-aje)e"'dxdy+y2f \du\^e'^~^dz 

for functions u G dom/ig = ^^(Mg) where ^^(Mg) is the closure of (^^(Mg) 
in the norm + hir{u)Y^'^. Here, \dyu\^ and \du\^ are evaluated in the {e- 

independent) metric of the exterior derivative of ^(a;, •) and u on T*F and T*Vk, 
respectively. 

The spectrum of ifg is again purely discrete (since is compact) and will be 
denoted by Xk{s), written in ascending order and repeated according to multi- 
plicity. By the the min-max principle we have 

Afc(e)=inf sup (2.5) 

«GLfe\{0} \\u\\e 

where the infimum is taken over all /c-dimensional subspaces of dom/i^, 
cf. e.g. [Di96|. 



We denote by the manifold F with metric e'^dy'^ and the first Dirichlet 
eigenvalue of F by Ai = Af (F) > 0. Since the lowest eigenvalue of F^ is Xi/e"^, 
we need a rescaling of the operator if^ in order to expect convergence to an 
^-independent limit operator. Therefore we set 

■■=H,-^ (2.6) 

and denote by the associated quadratic form. 
We first note that the operator is positive: 



Lemma 2.1. Suppose the smallness condition (ll.2j) is fulfilled, then > 0. 
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Proof. For u G dom = dom we have 



2 



+ ^{^dyu[x, ■)fp - Ai||m(x, OIID] (1 - rfx 



Applying the min-max principle for the first eigenvalue of the manifold F and 
Vfc, respectively, we conclude 

\\dyu{x, > Ai||n(x, and \\du\\l^ > X^"" {Vk)\\u\\l^. (2.7) 

Note that u ly^ lies in the quadratic form domain of A^^^. Using Assumption ()1.2p 
we see that qei^u) > 0. □ 

We set 

■■= WuWi + Qeiu) = \\u\\l + {h,{u) - 

Let us now formulate a simple consequence of the min-max principle ()2.5p which 
will be crucial in order to compare eigenvalues of operators acting in different 
Hilbert spaces (for a proof, see e.g. |EPn4[ Lemma 2.1]. Suppose that H, H' are 
two separable Hilbert spaces with the norms ||-|| and ||-||'. We need to compare 
eigenvalues and A^ of self-adjoint operators Q and Q' where Q > —A for some 
constant A > 0, with purely discrete spectra defined via quadratic forms q and 
q' onV cn and V C H'. We set \\u\\l := (1 + A)\\uf + q{u). 

Lemma 2.2. Suppose that J: T) — > T>' is a linear map such that there exist 
constants 61,62 > with 61 < 1/(1 + A + Afc) and 

\\uf < \\Ju\f + 6i\\u\\l (2.8) 
q{u) > q'{Ju) - 62\\u\\l (2.9) 

for all u eT). Then 

Afc > A^ - r]k 

where r]k is a positive function given by 

a A (Afc^i + ^2)(1 + A + Afc) 

Vk = v{Xk, 61, 62) ■■= — ...... X. — • (2.10) 

1 - (1 A Ak)6i 

In particular, ?7fc — >■ as 61,62 ^ 0. 
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3. Convergence of the eigenvalues: small vertex neighbourhoods 

In this section we consider a graph with straight edges approximated by an 
open set as defined in the previous section (the case of curved edges will be 
treated in the next section). We apply the abstract comparison result Lemma 
to our concrete problem in order to show an upper and a lower bound on \k{Qe) = 

Afe(ADj-Ai/£2 = Afc(£)-Ai/£2. 

Upper bound. We define the linear map Jo : dom go — ^ dom transmitting 
(eigen-)functions on the graph to functions on by 

(Jo/)(.):=e-/4^^"^^^'^' " = (^'^)^^^- (3.1) 

[0, z eVk 

where cp is the first normalised Dirichlet eigenfunction on the transversal direction 
F, i.e., 

Note that f\gj^ vanishes and therefore Jo/ ^ domg^ = ii^{M^). We begin with 
the verification of ()2.8|1 and ()2.9|1 . We have 

\\f\\l-\\ml = eY.aJ \f{x)\Ux = 0{e)\\f\\l (3.2) 

since ||(/?||_f = 1. Furthermore, 
qe{M) - qoif) 

= Y.\t^^ [ \ffdx + -J [ {\d,^\'-X,\v\')dy\f\\l-a,e)dx 

(3.3) 

Since (p is the eigenfunction with eigenvalue Ai the latter integral vanishes and 
therefore 

qeW) - qoif) = J^T^^Wf'Wl = 0{e)q,{f). (3.4) 

Applying Lemma f2.2l with A = we obtain 

Afc(£)-^< Afc(0) + O(£). (3.5) 

Lower hound. For the lower bound we have to work a little bit harder. We define 
Je : dom — > dom go by 

(J,m)j(x) := e'^'^{Nu{x) - p{x)Nu{x^)) (3.6) 

where 

Nu{x) := {u{x,-),ip) = I u{x,y)Tp{y)dy (3.7) 



10 OLAF POST 

is the expectation value of u{x, ■) G L2{F) corresponding to the lowest transversal 
eigenfunction ip. Here, x° depends on x and denotes the left resp. right endpoint 
of Ij if X is in the left resp. right half of Ij. Furthermore, p is a smooth function 
with < p{x) < 1, p{x) = near the mid point of Ij and p{x) = 1 near the 
boundary of Ij. Abusing the notation a little bit, x° also represents an element 



of dij. Since Jeu{x ) = 0, we have JsU G domgo- 



Again, we begin with the verification of ()2.8|) . First, we show the following 
estimate on higher transversal modes. 

Lemma 3.1. We have 

\\vr-\{v,^)\'<-^{\\dvr-x,\\vr) 

forv G Ti.^{F) where Aj are the Dirichlet eigenvalues of F. 

Proof. Since v — {v, ip)(p is the projection onto ip^, the min-max principle implies 

A2 

= UWdvf - Xi\{v,v)n = Uwdvf - A,||.f ) + ^{\\vf - \{v,^)n. 

A2 A2 A2 

Since F is connected, A1/A2 < 1 and we can bring the last difference on the LHS, 
divide by (1 — A1/A2) and obtain the desired estimate. □ 

The next lemma shows that under our main assumption, eigenf unctions do not 
concentrate at the vertex neighbourhoods: 



Lemma 3.2. Assume ()1.2|) then 



\u\\l . < 



2 r \ 



- Ar(Vfc) - Ai 
for all u G fi\M,) n T^HK.fc)- 

Proof. Using the second estimate in ()2.7|1 and the scaling of the metric fl2.3|l we 
have 



\u 



2 ^2 

2 



Ai 



2 



By our main assumption p.2|l . Ai/Af^(Vfe) < 1 and the result follows as before. 

□ 

Finally, we need the following lemma. 
Lemma 3.3. We have 

e-|Ar^x°)p < 0(e) OlHIi, - ^IkllU 
for all u G i-i^^Mi,) fl HMV^^k) o,nd x° G dIj, if Cj emanates from v^. 
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Proof. A standard Sobolev embedding theorem gives 

\Nuix')\' < \u{x\y)\''dy < c,[\\dury^ + ||«||y 

for some constant Ci > (note that F = djVk). Now by the scahng of the metric 
on Vfc 



and the result follows from the preceeding lemma. 
Now, we want to consider the norm difference 



IK 



□ 



W Z W J w z 
Me - \\-JeU\\o 



ieJ 



\u\ 



\Nu{x) - p{x)Nu{x'^)\'e"'dx 



The first sum can be estimated by 0{e'^)qs{u) using Lemma f3.2[ For the second, 
we use 



{a + bf > {l-6)a^--b^, 





6>0 



(3i 



and obtain as upper bound 



\u\ 



-(1-5) / \Nu{x)Ye'^dx + 



\p{x)\^dx max \Nu{x 



0m2 



< 



|m(x, Of - |(m(x, I 



e'^dx 



+ 



5 



1 — ttjE 



a^e \\u 



0\|2 



using Cauchy-Schwarz. Applying Lemma the scaling of the metric on F 
in ()2.3j) . Lemma [3.31 and setting 5 = e^^"^, we end up with the estimate 



\u\ 



\JM\l<Ois'/')\\u\\l,. 



(3.9) 



For the quadratic form difference we have 
qoiJeu) - qe{u) 



< 



y^e-^\\N{d^u)-p'Nu{x'^)\\l-—^ I WdM^rWe'^dx 

^ — ^ 1 — I J . 



ieJ ^ 
where the terms of order has been estimated with ()2.7|) . Using 

[a^hf <{\^S)o? ^h?, 0<(5<1, 





(3.10) 
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with 5 = e^^"^, Cauchy-Schwarz for \N{dxu{x, ■)p < \\dxu{x, and Lemma IH7^ 
we obtain 

qoiJeu) - Qeiu) < Oie'/')\\u\\l,. (3.11) 
Applying Lemma \'2.'2\ again (with A = 0), we obtain 

XkiQe) = Afc(5) - ^ > A,(0) - 7]k- (3.12) 

Here, r]k = 0(e^/^) using ()3.9|1 . (jH.llj) and the upper estimate Xk{Qe) < Afc(O) + 
0{e) = 0(1) from 

4. Curved edges 

Let us now consider a curved quantum wave guide embedded in (more 
general embeddi ngs can be treated similarly). Such spaces have already been 



analysed e.g. in ES89 or |DE95j . We only consider a single edge here since 



one can easily replace the edge estimates without curvature by the appropriate 
estimates with curvature in the previous section^ (cf. Remark 14.11 for the precise 
assumptions on the curvature). The convergence of the discrete spectrum of an 
infinite curved quantum wave guide has already been established in |DE95j using 
perturbation arguments and an asymptotic expansion (cf. also |KaP88j where the 
asymptotic of the first Dirichlet eigenvalue of a ^-neighbourhood of a finite length 
curve in was treated). Here, in contrast, we use the variational arguments of 
Lemma [2.21 which are somehow simpler (the price being a weaker result). 

Definition of the approximating space. Suppose that 7: / — > is a smooth 
curve (e.g. is enough) with bounded derivatives parametrised by arclength 
(i.e. the tangent vector 7(x) has unit length for all x G /). Suppose that either 
7 is a closed curve (/ = S^) or has two ends (/ = (0, 1)). 

We introduce the ^-neighbourhood of the curve given as the image of the 
parametrisation 

{x,y) I — > 'y{x) + ey n{x) ^ ' ' 

where n(x) := (72(2^), —71 (a;)) is orthogonal to the tangent vector 7(0;) and F = 
(— 1, 1). Define the signed curvature by 

n ■= 7i72 - 727i (4.2) 

and suppose that < e < Eq := 1/(2||k||oo) where ||k||oo denotes the supremum 
of |k(x)|, X G /. We assume in addition that \1/ is a diffeomorphism. 



""^More precisely: one has to show the estimates of this section for the metric 

(1 - ea)^(l + ey n{x)Ydx^ + e^dy'^ 

instead of the metric defined in 1)4. 3|) in order to take into account the shortened edges due to 
the embedding. To keep the notation manageable we omit this fact here. 
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Remark 4.1. Suppose we consider an embedded graph Mq with curved edges 
Cj with curvature kj. Besides the assumption that the parametrisation (j4.H) 
is a diffeomorphism for each edge cj we need the additional hypothesis that 
supp Kj is contained in the open interval /-,-, i.e., that the curvature vanishes in a 
small neighbourhood of the adjacent vertices. Otherwise one needs to modify the 
scaling property of the vertex neighbourhoods Ve^k'- they cannot be homothetic 
to a fixed set Vk if the edge is curved up to the vertex Vk- 

Denote by the pull-back of the Euclidean metric via i.e., Qg := \l/*5'euci- 
A straightforward calculation shows that 

= {1 + ey K{x)ydx'^ + e'^dy'^. (4.3) 

We denote by the manifold I x F with metric and by the same manifold 
with the product metric 

Ijs = dx^ + e^dy"^. 

For computational reasons, it is much easier to deal with the latter metric so we 
introduce the unitary transformation 

u I — > {I + ey K{x)f''^u. ^ ■ ' 

1 /2 

Note that det Qe = e{l + ey k{x)) > is the density of the metric and 

~l/2 

deiQe = £■ For the rest of the section, we will work in the Hilbert space 

Definition of the operator on the thickened set. We want to consider the Dirichlet 
Laplacian on f/^. Its quadratic form is ||(i-u||^^, u G iiMUg) (we could also allow 
other boundary conditions at dl x F). The transformed quadratic form is given 
by 

K{u) := \\d^*u\\l = Mil + ey k)-'/'u)\\1, u E H^U,) = n\U,) 



A straightforward calculation already performed at various places (e. g. |ES89 
or jDE95]) yields 

he{u) = I I . —\d.j,u\^^\\dyu\^^Kg\u\^edydx (4.5) 

where the curvature induced potential is given by 

^^""'^^ A{1 + ey ^ 2{1 + ey A{1 + ey k)^' ^ ' 

Note that Kir, < e < eo, is bounded from below by a constant — A^^, Ag^ > 
depending only on the supremum of and eQ. Using the first estimate 

in ()2.7|) we see that 

qeiu) ■.= he{u) -^\\u\\}j^ (4.7) 



14 OLAF POST 

is bounded from below by — A£,J|n||^ . Tfierefore, 

■■= Qeiu) + (A,, + l)\\u\\l^ 
defines a norm on tlie quadratic form domain ii}{Ue). 

Definition of the limit space and operator. Finally, we define tlie limit operator 
Qq. Clearly, Qq will act in the limit space Tio := L2{I). As usual, we define Qo 
via its quadratic form 

,2 



^o(/):= J^[\fr-^\f\'\dx. (4.8) 

Again, 

ll/l&i :=go(/)+(^ + l) 11/11? 

defines a norm on the quadratic form domain 'fi^[I). Note that K^{x,y) = 
-K{xY/4: + 0{e) as £ ^ 0. 

Spectral convergence. We want to show the following spectral convergence. From 
its proof it is straightforward to show Theorem ll.ll in the general case of branched 
quantum wave guides with curved edges. 

Theorem 4.2. Denote by Xk{£) the k-th Dirichlet eigenvalue of the curved quan- 
tum wave guide Us. Then 

Afc(£)-^ = A,(0) + O(£), e^O, 

where Afc(O) denotes the k-th eigenvalue of Qq = —d^/d^x — /A. 

Here, Ai = 71^/4 is the first Dirichlet eigenvalue of F = (—1, 1). As before, we 
establish an upper and a lower bound on \k{s)- 

Upper hound. We define the transition operator Jq as in (jH.lll on the edges (here, 
m = 1). Clearly, we have 

11/11^=11^0/11.^ 
since ip is supposed to be normalised. In addition. 



e dy dx 



II [( (TT^ - 1) i/vi^ + i/r^^(i.f - A. I.P) + {K. . ^) l/.P 

which can be estimated by C'(£:)||/||oi where 0{e) depends only on k (and its 
derivatives) (remember that ip is the first Dirichlet eigenfunction on F with 
eigenvalue Ai). Applying Lemma 12.21 yields the desired upper estimate with 

Vk{e) = 0{e). 
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Lower bound. The lower bound is again a little bit more difficult. We define the 
transition operator by 



{Jeu)ix) ■.= e^^^Nu{x) 



(4.9) 



where Nu is the transversal expectation value of u with respect to if, cf. Equa- 
tion ()3.7|) . Applying Lemma ITT] for v = u{x, ■) we obtain the estimate 



\u\ 



\JsU\\o < 



A2 — Ai J J 



— (I^y-ul^ — Ailwp) edydx 



< 



A2 — Ai 

The quadratic form estimate is given by 



qo{Jeu) - qe{u) 



{\{d^u{x,-),^)\ 



1 



l + eyK) 



Aikr)+^(|txr-|Kx,-),¥^)P) 



\u\ 



K 



e dy dx. 



As before, we estimate the ffist difference using Cauchy-Schwarz. The second 
difference is negative (cf. ()2.7|) ). The third difference is small due to Lemma f3. II 
The forth difference is also small since = — + 0{e). Therefore, we end 
up with an upper estimate given by (^(e:)!!^!!^ ^. Applying Lemma f2. 21 once more, 
we obtain the desired lower estimate on Xkie)- Using also the upper estimate we 
see that rik{e) = 0{e). 



5. Examples 

In this section, we want to comment on the smallness condition fll.2|) and give 
examples where this condition holds or fails. To simplify the presentation, we 
assume that Mq C M^. 

First, we show, that the condition can always be fulffiled, provided the vertex 
neighbourhood is small enough. Suppose that we start with the 1-neighbourhood 
denoted by Vfc(O), i.e., we set e = 1 and regard the unsealed vertex neighbour- 
hood Vfc. Remember that we have assumed that the curvature vanishes near 
the vertices, therefore Vk{0) is bounded by straight lines. Then we deform Vfc(O) 
smoothly in order to obtain a family Vk{T), r > 0, shrinking to the graph, but fix- 
ing the boundary parts djVk{T) = djVk{0), j G Jk, where the edge neighbourhoods 
touch (cf. Figure El). As in |P03l Sec. 7] we can show that the ffist eigenvalue of 
the Laplacian on Vfc(r) with Dirichlet boundary conditions except on the fixed 
boundary parts 9j\4(r), where we impose Neumann boundary conditions, tends 
to 00, i.e., 

Af^(Vfc(r)) ^00 as r ^ 00. 



16 OLAF POST 




Figure 3. The original vertex neighbourhood Vfc(O) (hght grey) 
and the shrunken vertex neighbourhood Vfc(r) (dark grey). 

Therefore there always exists a fixed r G (0, oo) such that Vk := Vfc(r) satis- 
fies (II .21) . Fixing this shrinking parameter r, we proceed with the definition of 
Me as in Section |21 

An example not satisfying the smallness assumption. Let us briefiy give an ex- 
ample of a vertex neighbourhood not satisfying (|1.2j) . For suitable vertex neigh- 
bourhoods (e.g. arising from the ^-neighbourhood of a graph) we will show the 
existence of an eigenvalue below the threshold Xi/e^ = tt^/ (4e^) (cf. also |SRW89j 
and |ABGM9T] for the case of an e-neighbourhood of a vertex with four infinite 
edges emanating (a "cross"); in the former reference one can also find a contour 
plot of the first eigenfunction). Therefore, the conclusion of Theorem II. II is false, 
i.e., dig) fails. 

The existence of such an eigenvalue below the threshold can easily be estab- 
lished by inserting an appropriate trial function in the Rayleigh quotient. We 
consider a graph with one vertex and three adjacent edges of length i and denote 
its ^-neighbourhood by M^. We decompose into three rectangles f/^j- and 
three sets A^j as in Figure El 

On the rectangle Usj we use the coordinates < x < £ and —e<y<e 
where x = corresponds to the common boundary with A^j. We extend these 
coordinates from Ue,j onto j and define 

u{x,y) ■=e~^l'^x{x) cosily) 

as a test function on each of the three sets Usj U A^j. Here, x{^) = 1 for x < 
(i.e., on A^j), x(x) = cos(7rx/(2£/t)) for < x < ne and x{^) = for e < x < i 
where k > is some constant to be specified later. Although u is not differentiable 
across the different borders (but continuous), it still lies in the quadratic form 
domain n^{M,). 
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Figure 4. A simple trial function supported in a neighbourhood 
of the vertex has an eigenvalue below the threshold Xi/e"^ = 
7rV(4£2). 



A straightforward calculation yields 

IMi^IIm^ vr^ / Skcosq; + Svr^ sina — 16fi; 



\u 



12 

I Me 
|2 



4^2 



((Stt^ — 4) cos a + Svr^K sin a + 8)hi/ Ae^ 



(5.1) 



This quantity is negative for all < a < O.QSvr if we choose e.g. k = 3. In 
particular, there exists a negative eigenvalue of Aj^^^ — Xi/e^ of order e"^, and 
Condition (|1.2j) fails here for any choice of vertex neighbourhoods V^^k since the 
conclusion of Theorem II. II is false. Note that the vertex neighbourhoods V^.u are 
not uniquely determined. One could enlarge V^^k at each edge emanating by a 
cylinder of length ae taken away from the corresponding edge neighbourhood. 

If — \i/e^ has negative eigenvalues it is not clear whether its appropriately 
scaled eigenvalues converge to eigenvalues of an operator on the graph Mq. The 
dependence of the leading order on the angle a in ()5.1|) indicates that the limit 
should depend on the angles of the edges meeting at a vertex. 
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BRANCHED QUANTUM WAVE GUIDES WITH DIRICHLET 
BOUNDARY CONDITIONS: THE DECOUPLING CASE 



OLAF POST 

Abstract. We consider a family of open sets which shrinks with respect to 
an appropriate parameter e to a graph. Under the additional assumption that 
the vertex neighbourhoods are small we show that the appropriately shifted 
Dirichlet spectrum of converges to the spectrum of the (differential) Lapla^ 
cian on the graph with Dirichlet boundary conditions at the vertices, i.e.. a 
graph operator without coupling between different edges. The smallness is ex- 
pressed by a lower bound on the first eigenvalue of a mixed eigenvalue problem 
on the vertex neighbourhood. The lower bound is given by the first transversal 
mode of the edge neighbourhood. We also allow curved edges and show that 
all bounded eigenvalues converge to the spectrum of a Laplacian acting on the 
edge with an additional potential coming from the curvature. 



1. Introduction 

Graph models of quantum systems can often be used to describe in a simple 
way some important aspects of the behaviour of a quantum system. Although 
such models are simple enough to be solvable (because they are essentially 1- 
diniensional) they still have enough structure to model real systems. Riiedenberg 
and Scherr [RuS53] used this idea to calculate spectra of aromatic carbohydrate 
molecules. Nowadays the rapid technical progress allows to fabricate structures 
of electronic devices where quantum effects play a dominant role. Graph models 
like quantum graphs (also called metric graphs) can often l)e viewed as a good 
approximation of such structures. From the mathematical point of view these 
models were analysed first thoroughly in [ES89] , for recent developments, bibli- 
ography and further appUcations see [DE95], [KoS99], [Ku02] or [Ku04]; note that 
[KaP88] also calculated the eigenvalue asymptotic of a tubular e-neighbourhood 
of a curve. 

A quantum (or metric) graph is a graph where we associate a length to each 
edge. A natural operator acting on such graphs is given by a self-adjoint extension 
of —tP/dx'^ on each edge. We will call such a self-adjoint extension a Laplacian on 
the ((luantum) graph. Note that the Laplacian on a discrete graph is a difference 
operator on f-2(K) rather than a differential operator acting in L2{ej). Here, 
K labels the vertices and ./ the edges Cj, j 6 J, of the graph. A detailed overview 
on this wide field can be found in [Ku04] or [KoS99]. 
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A natural question is in wliat mathematical sense a quantum graph Mq can be 
approximated by a more smooth space M^. One is interested what Laplacians 
on Mg occur as limit operators from operators on M^. More significantly, Mj 
could be the £-neighbourhood of an embedded graph Mq c K" or a manifold 
shrinking to Mq as e ^ 0. We call such approximating spaces branched quantum 
wave guides. Recently, spectral convergence in the case of a bounded open set 
AI^ with Neumann boundary condition has been established in [RScOl], [KuZOl] 
and [KuZ03]; for an approximation by manifolds see [EP04]. All these examples 
have in common, that the lowest eigenmode of the transversal direction is with 
constant eigenfunction. In this case, the limit operator is the Laplacian on the 
graph with Kirchhoff boundary condUions, i.e., a function / in the domain of the 
Kirchhoff Laplacian is continuous at each vertex and satisfies 

E /iK) = 0, keK. (1.1) 

In addition, the spectral convergence holds independently of a given embedding 
of the graph. In particular, the convergence is independent of the ciu'vature of 
the embedded edges. 

The case of an approximation by Dirichlet Laplacians on an open set was 
first treated liouistically in [RuS53]. This case is harder to analyse since the first 
transversal eigenvalue equals Af (F^) = Xi/e^ (Ai > 0), i.e., it is of the order e^^ 
if £ denotes the radius of the cross section = (— e, e) of the approximating set 
Mj. A rescaling is necessary, and first order terms of the metric (cf. (4.3)) like 
the curvature become important. In particular, the curvature of the (embedded) 
edge enters in the Umit operator as an additional potential. 

Main result. Assume that Mq c is a finite graph. Our aim in this note is to 
show the spectral convergence of the Dirichlet Laplacian on an approximating 
open set D Mi-,. We suppose that can be decomposed into neighbour- 
hoods C4j- of the edges ej and into neighbourhoods ICjt of the vertices Vk of Mq 
(cf. Figure 1). We assume that Vs,k is e-homothetic to a fixed set Vk- The precise 
definition will be given in Section 2 and 4. Our basic assumption is that the 
vertex neighbourhoods Ve^k are small, i.e., that 

Ar(H) > A?(F) = Ai (1.2) 

where A™(Vt) is the lowest eigenvalue of the Laplacian A™ of 14 with Dirichlet 
boundary conditions on daVk (i.e., on the boundary induced from the original 
boundary of M^) and Neumann boundary conditions on djVk, j £ .A-, (i.e., on 
the parts where the adjacent edge neighbourhoods labeled by j 6 Jk emanate, 
cf. Figure 2). Furtluaniore, F = (—1, 1) and therefore Ai = 7r^/4. We comment 
on this condition in Section 5. 
Our main result is 
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Figure 1 . Decomposition of the graph neighbourhood AI^ (gi'f^y) 
of the graph i\/o into edge and vertex neighbourhoods and V^^k- 




Figure 2 . The scaled vertex neighbourhood T4 with the boundary 
part doVk coming from the original boundary and the boundary 
part djVk where the edge ej emanates. 

Theorem 1.1. Suppose that is an open neighbourhood of a finite graph 
Mq C satisfying the smallness assumption (1.2) on each vertex neighbour- 
hood. Denote by Xkie) the k-th eigenvalue of the Dirichlet- Laplacian > 
(counted with respect to multiplicity). Then 

A.(£) - ^ ^ A,(0). e^O (1.3) 
where Xk{0) denotes the k-th eigenvalue of 
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with A^r" — —(p/dxj being the Dirichlet Laplacian on the edge ej = Ij = 
and Kj being the curvature of the embedded edge C (cf. (4.2) j. 

Note that the smallness assumption at the junctions V^^i^ imphes that the limit 
operator decouples, i.e., the limit operator is the direct sum of operators acting 
on a single edge. In the case of the Neumann Laplacian on decoipliiig occurs 
if the area of the edge neighbourhood decays faster than tlie area of tlie vertex 
neighbourhood; e.g., if the latter scales in each direction of the order s"' with 
< Q < 1/2; the vertex neighbourhoods are large ol)stacles seen from the edge 
neighbourhoods (cf. [KuZ03] or [EP04]). In the case of Dirichlet l)oun(lary con- 
ditions, in contrast, decoujjling already occurs when the vertex neighbourhoods 
scale with e, i.e., even when tlie edge neighbomhood volume (which is of order 
e) decays slower than the vertex neighbourhood volume (of order e^). 

We also show in Section 5 that the usual £-neighbourhood := {z e 
I d{z, Mo) < £ } does not satisfy our hypothesis since the leading order of 
the lowest eigenvalue is at most /i/e^ with yu < Ai. Therefore, — Ai/e^ has 
a negative eigenvalue tending to — oo (cf. also Lemma 2.1) and the conclusion 
of Theorem 1.1 fails. In particular, there is no limit operator on the graph (us- 
ing the simple shift A — > A — Ai/e^), and the suggestion in [RuS53], that the 
limit operator is the Kirchhoff Laplacian on M^, is false (cf. also [Kti02, Sec. 2.1 
and 3.2]). Note that Ruedenberg and Scherr implicitly assumed that the lowest 
eigenfunctiou does not concentrate around the vertex which is the case as we will 
see in the last section. 

The spectral convergence of a single curved quantimi wave guide has already 
been shown in [DE95] and [KaP88] using perturbation methods. Our proof only 
uses variational methods and is a simple adaption of [EP04], [KuZOl, KuZ03] or 
[RScOl], where one compares Rayleigh quotients. 

The paper is structured as follows: In the next section we define the required 
spaces and operators in the case of straight edges. Section 3 is devoted to the 
proof of Theorem 1.1 in this case. In Section 4 we provide the necessary changes in 
order to prove the result with curved edges. Section 5 contains some explautation 
on the smallness condition (1.2), and examples where this condition holds or fails. 



2. Preliminaries 

In this section we define the limit space and the approximating space together 
with the associated operators. We hrst consider straight edges without curvature, 
i.e., Kj = 0. In Section 4 we also allow curved edges. 

Definition of the limit space. Let Mg be a finite connected graph with (metric) 
edges Cj, j & J {cj isometric to an open interval /_,■ — {O.f.j)) and vertices 'u^, 
k e K. We denote the set of aU j € J such that Cj emanates from the vertex Vk 



BRANCHED QUANTUM WAVE GUIDES 



TT 










i 


/ 


a \ V 










~"""-Ae.2 














Figure 4. A simple trial function supported in a neiglibomhood 
of the vertex has an eigenvalue below the threshold Xi/s'^ = 
7rV(4£2). 



A straightforward calculation yields 



\\du\\l 



IT 

4£2 



16k 



/ 8k cos a + 37r^ sin a ■ 
V {{3tt'^ — 4) cos a + 37r^K sin a + 8)k'' 4e^ 



(5.1) 



This quantity is negative for all < a < 0.937r if we choose e.g. k = 3. In 
particular, there exists a negative eigenvalue of A]^,^ — Ai/e^ of order e"^, and 
Condition (1.2) fails here for any choice of vertex neighbourhoods l^jt since the 
conclusion of Theorem 1.1 is false. Note that the vertex neigliboiulioods Vs,k a-re 
not imiquely determined. One could enlarge V, ^ at each edge emanating by a 
cylinder of length as taken away from the corresponding edge neighbourhood. 

If A^^ — Ai/e^ has negative eigenvalues it is not clear whether its appropriately 
scaled eigenvalues converge to eigenvalues of an operator on the graph Mq. The 
dependence of the leading order on the angle a in (5.1) indicates that the limit 
should depend on the angles of the edges meeting at a vertex. 
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Figure 3. The original ^xTtcx iK!iglil)ourhoo<i Vk{0) (light grey) 
and the shrunken vertex neighbourhood 14 (r) (dark grey). 



Therefore there always exists a fixed r £ (0,oo) such that Vk := Vjfe(r) satis- 
fies (f.2). Fixing this shrinking parameter r, we proceed with the definition of 
Mj as in Section 2. 

An example not satisfying the smallness assumption. Let us briefly give an ex- 
ample of a vertex neighbourhood not satisfying (1.2). For suitable vertex neigh- 
bourhoods (e.g. arising from the e-neighbourliood of a graph) we will show the 
existence of an eigenvalue below the threshold Ai/e^ = 7r^/(4e^) (cf. also [SRW89] 
and [ABGM91] for the case of an e-neighbourhood of a vertex with four infinite 
edges emanating (a "cross" ) ; in the former reference one can also find a contour 
plot of the first eigenfunction) . Therefore, the conclusion of Theorem 1.1 is false, 
i.e., (1.2) fails. 

The existence of such an eigenvalue below the threshold can easily be estab- 
hshed by inserting an appropriate trial function in the Rayleigh quotient. We 
consider a graph with one vertex and three adjacent edges of length t and denote 
its £-neighbourhood by M^. We decompose Mj into three rectangles U^j and 
three sets A^j as in Figure 4. 

On the rectangle Uej we use the coordinates Q < x < t and —e<y<£ 
where s = corresponds to the common boundary with A^j. We extend these 
coordinates from Ue,j onto A^j and define 

u{x, y) := s-^^\{x) cos(^</) 

as a test function on each of the three sets U^j U A^ j. ffere, xi^) = 1 for a: < 
(i.e., on A,j), x(,-':) = cos(7r,T/(2cK)) for < x < ne and x(-'') = for e < < ^ 
where k > is some constant to be specified later. Although u is not differentiable 
across the different borders (but continuous), it still lies in the quadratic form 
domain fO^iMe). 
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by Jfc. The Hilbert space associated to such a graph is 
H:=L2(Mo) = 0L2(/,) 

which consists of all functions / with finite norm 

ii/ii§ = ii/iiL„ = Eii/^iil = E / \m\"dx. 

Definition of the limit operator. We define the limit operator Qo via the quadratic 
form 

*{/) :=Eii/iiil = E / if'M'dx 

for functions / e C~(Mo) = ^jC^ilj) (with compact support). The form 
closure of go (also denoted by go) is the extension of go to the closure of the space 
of all such functions in the norm 

II/II0.1 ll/llo + *(/) 

(see [K66, Chapter VI], [RS80] or [Da96] for details on quadratic forms). Note 
that 

d„m,^„ = 07^i(/,). 

Remember that fl^I) is the closure of C^{I) w.r.t. the norm {\\f\\^ + Wf'fV^^- 
The associated self-adjoint, non-negative operator Qo is given by 

Qo = 0A° (2.1) 

where A^. > denotes the self-adjoint operator —d?/dx^ on /, with Dirichlet 
boundary conditions. The spectnmi of Qo is purely discrete and will be denoted 
by Ai,(0), written in ascending order and repeated according to nniltiplicity. 

Definition of the approximating space. We now describe the family of open sets 
{M^)g, < e < So, approximating the graph Mo as e ^ 0. For convenience 
only, suppose that Mo is embedded in (au abstract definition of in the 
general case will be given soon). Assume that we can decompose A/^ into open 
sets Us.j containing those points x 6 ej with d{x,de.j) > ci_,e/2 for some real 
number Oj < l/co and K ^ 3 'I'k such that the union of their closures equals 
M^. Here, the edge neighbourhood U^j is isometric to I^j x (both equipped 
with the Euclidean metric) where /jj := (0, (1 — aje)ij) and = {—e.e) is the 
scaled cross section. Furthermore, we assume that the vertex neighbourhood V^^j, 
is £-homothetic to a fixed open set Vk. Using a simple coordinate transform we 
have therefore the isometrics 



{Ue,j,geuc\} ^ {I3 X F,ge) and {Ve,k,geuc\}^ {Vk,ge) (2.2) 
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where 

ffj = (1 - eajYdx^ + e'^dy'^ and = e^g (2.3) 

are the metrics on the edge resp. vertex neighbourhood. Here, F = (— 1, 1) and 
g is the Euchdean metric on 14 . In the sequel we use this change of coordinate 
transform without mentioning. Note that the slightly shortened edge neighbour- 
hood is necessary in order to have an embedding for the edge and the vertex 

neiglibourhood. 

Althougli we are mainly interested in the embedded situation as described 
above, we prefer the following abstract setting in order to keep the notation 
of [EP04] and recognise the important geometric objects (not depending on any 
embedding). For each j e J we let U^j be the Riemannian manifold [Ij x F, g^) 
where is given as in (2.3). Here F is the interior of a compact, connected m- 
dimensional manifold (m > 1) with metric denoted by dy'^ having purely discrete 
Dirichlet spectrum with first eigenvalue Ai > 0. 

Furthermore, we denote by Ve^k the Riemannian manifold (T4, g^) with g^ = e^g 
where g is a metric on 14 . We assume that 

= doVk U U djV^, (2.4) 

i.e., the boundary of Vk has as many boundary parts djVk isometric to F as 
edges emanate from Vk and doVk is the closure of dVk \ \Jj djVk (cf. Figure 2). 
Furthermore, we assume that the metric on Vk has product structure g = dx^+dy^ 

near djVk. 

We can define an abstract manifold by identifying the api)ropriate boimdary 
parts according to the graph Mo. Note that a smooth structure on and also 
a smooth metric g^ of the form (2.3) in the respective charts exist since is 
diffeomorphic to a product (0, 1) x F in a neighbourhood of each djVk on both 
sides of djVic, i.e., on 14 and Uj. Strictly speaking we should introduce another 
chart for each j e Jk and k e K covering djVk in order to define the smooth 
structure properly. But since we only use integrals over M^, a cover up to sets 
of measure is enough. The resulting manifold has dimension d = m + 1. Note 
that Me need not to be embedded in any space, but the embedded case described 
above is also covered by this setting. 

The associated Hilbert space is 
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Lower bound. The lower bound is again a little bit more difficult. We define the 
transition operator by 

{Jeu)(x) ■.= e^l'^Nu{x) (4.9) 

where Nu is the transversal expectation value of u with respect to ip, cf. Equa- 
tion (3.7). Applying Lemma 3.i for v = u{x, •) we obtain the estimate 

A2 — Ai 

The quadratic form estimate is given by 

- Ail^p) + '^{\u\' - \(u{x,-),v)\') - \u\\'^ + K,)]edydx. 

As before, we estimate the first difference using Cauchy-Schwarz. The second 
difference is negative (cf. (2.7)). The third difference is small due to Lemma 3.1. 
The forth difference is also small since = — + 0(e). Therefore, we end 
up with an upper estimate given by 0(e)|ju|lj j. Applying Lemma 2.2 once more, 
we obtain the desired lower estimate on Xkie:)- Using also the upper estimate we 
see that ^(e) = 0(e). 

5. Examples 

In this section, we want to comment on the smallness condition (1.2) and give 
examples where this condition holds or fails. To simplify the presentation, we 
assume that Mq C M^. 

First, we show, that the condition can always be fulfilled, provided the vertex 
neighbourhood is small enough. Suppose that we start with the 1-neighboiu'hood 
denoted by 14(0), i.e., we set e = 1 and regard the unsealed vertex neighbour- 
hood 14. Remember that we have assiuned that the ciu'vature vanishes near 
the vertices, therefore 14(0) is bounded by straight lines. Then we deform 14(0) 
smoothly in order to obtain a family Vt(r), r > 0, shrinking to the grajjli. but fix- 
ing the boundary parts djVk^r) = &,14(0). j G Jk, where the edge neighbourhoods 
toiicli (cf. Figure 3). As in [P03, Sec. 7] we can show that the first eigenvalue of 
the Laplacian on 14 (t) with Dirichlet boundary conditions except on the fixed 
boundary parts 9,14(r), where we impose Neumann boundary conditions, tends 
to 00, i.e., 

Ai'^'(14(t)) — > 00 as r ^ 00. 
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is bounded from below by — Ae„||«||?, . Therefore, 

\\u\\i, -q,{u) + iA,„ + l)\\u\\l 
defines a norm on tlie cjuadratic form domain W^(C4). 

Definition of the limit space and operator. Finally, we define the limit operator 
Qo- Clearly, Qo will act in the limit space Ho := L2{I). As usual, we define Qo 
via its quadratic form 

?o(/):=^[l/f -^l/P]rf^. (4.8) 

Again, 

:=*(/) + (^ + 1)11/11? 

defines a norm on the quadratic form domain n}[I). Note that Kc{x,y) — 
-«;(a;)74 + 0(e) as £ ^ 0. 

Spectral convergence. We want to show the following spectral convergence. From 
its proof it is straightforward to show Theorem 1.1 in the general case of branched 
quantum wave guides with curved edges. 

Theorem 4.2. Denote by \k{£) the k-th Dirichlet eigenvalue of the curved quan- 
tum wave guide Ue. Then 

Xk{e)-^ = Xk{0) + O{e), e ^ 0, 

where Xk{0) denotes the k-th eigenvalue of Qo = —iP/(fx — 

Here, Ai = 7r^/4 is the first Dirichlet eigenvalue of F = (—1, 1). As before, we 
establish an upper and a lower bound on Aji,(£). 

Upper bound. We define the transition operator Jq as in (3.1) on the edges (here, 
m= I). Clearly, we have 

ll/llo = lko/||? 
since ip is supposed to be normalised. In addition, 

QeW) - qo(f) = 

which can be estimated by O(e)||/||o i where 0{s) depends only on k (and its 
derivatives) (remember that ip is the first Dirichlet eigenfunction on F with 
eigenvalue Ai). Applying Lemma 2.2 yields the desired upper estimate with 
Vkis) = 0{e). 
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which consists of all functions u with finite norm 



f \u\^(l-aje)£"'dxdy-\-Y f \u\'e'^dz 



where dy and dz represent the natural measures on F and 14, respectively. 

Definition of the operator on the manifold,. The operator on the thickened space 
we are considering will be the Dirichlet Laplacian on M^, i.e., = > 0. 
The corresponding quadratic form he is given by 



M«) = 5^||(to||?,^ + ^|IHI|,^ 

= E/ [ n ^ .2 \dM' + -2M]{l-a,e)e^dxdy + Y, ( \du 



for functions u e dom/ij = fi^{Me) where 7l!^(Mj) is the closure of C^{Ms) 
in the norm (||«||^ + he{u)y^^. Here, {dyul"^ and \du\'^ are evaluated in the (e- 
independent) metric of the exterior derivative of u{x, ■) and u on T*F and T*14, 

respectively. 

The spectrum of is again purely discrete (since is compact) and will be 
denoted by Xk{s), written in ascending order and repeated according to multi- 
plicity. By the the min-max principle we have 

Afc(£) = inf sup (2.5) 

^1= u€Lt\m \m\e 

where the infimum is taken over all A)-dimensional subspaces Lj, of dom/i^, 
cf. e.g. [Da96]. 

We denote by F^ the manifold F with metric e^dy'^ and the first Dirichlet 
eigenvalue of F by Aj — Af (F) > 0. Since the lowest eigenvalue of F^ is \i/e^, 
we need a rescaling of the operator Hg in order to expect convergence to an 
£-independent limit operator. Therefore we set 

Qe ■.= He-^ (2.6) 

and denote by Qe the associated (quadratic form. 
We first note that the operator Q, is j)ositive: 

Lemma 2.1. Suppose the .smallness condition (1.2) is fulfilled, then > 0. 
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Proof. For u € doin = dom we have 



+ \{\\dyu{x, -Wi, - Xi\\u{x, Olll)] (1 - «,£)£'" d-r 



Applying the min-max principle for the first eigenvalue of the manifold F and 

14, respectively, we conclude 

\\dyU(x,-)\\%>Xi\Hx,-)\\l and >A™(y,)||«||^^. (2.7) 

Note that u\y^. lies in the quadratic form domain of A™. Using Assumption (1.2) 
we see that qe{u) > 0. □ 



We set 



Ml ■■= \\u\\l + qs{u) = \\u\\l + [Kin] - ^\\u\\l) 



Let us now formulate a simple consequence of the miu-max jjrincijjle (2.5) which 
will be crucial in order to compare eigenvalues of operators acting in different 
Hilbert spaces (for a proof, see e.g. [EP04, Lemma 2.1]. Suppose that Ti., H.' are 
two separable Hilbert spaces with the norms ||-|| and We need to compare 
eigenvalues Aa, and A'^, of self-adjoint operators Q and Q' where Q > —A for some 
constant A > 0, with purely discrete spectra defined via quadratic forms q and 
q' onVcn and V C H'. We set \\u\\l := (1 + A)\\u\\^ + q{u). 

Lemma 2.2. Suppose that J: V — > V is a linear map such that there exist 
constants 61,62 >0 with 61 < 1/(1 + A + A^) and 

\\uf < \\Ju\f + 5i\\u\\l (2.8) 
q{u)> q'(Ju)-S2\\u\{i (2.9) 

for alluEV. Then 

Afc > A^ - % 
where rjn is a positive function given by 

„ r.(\ fi (A,gl+^2)(1+A + A,) 



In particular, % — > as 5i, 52 — > 0. 
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Remark 4.1. Suppose we consider an embedded graph Mq with curved edges 
ej with curvature Kj. Besides the assumption that the parametrisation (4.1) 
is a diffeomorphism for each edge Cj we need the additional hypothesis that 
suppKj is contained in the open interval Ij, i.e., that the curvature vanishes in a 
small neighbourhood of the adjacent vertices. Otherwise one needs to modify the 
scaling property of the vertex neighbourhoods Vs^k- they carmot be homothetic 
to a fixed set Vj, if the edge is curved up to the vertex Vk- 

Denote by the pull-back of the EucUdean metric via 5', i.e., g^ := 5'*(?euci- 
A straightforward calculation shows that 

g, = {l + ey K{x)fdx'^ + e^dy^. (4.3) 

We denote by the manifold I x F with metric g^ and by Ug the same manifold 
with the product metric 

Ije = dx^ + e^dy^. 

For computational reasons, it is much easier to deal with the latter metric so we 
introduce the unitary transformation 

u ^ (1 + ey K{x)Yriu. ^ ' 

Note that detge'^ = e(l + ey k{x)) > is the density of the metric g^ and 
detge'^ — e. For the rest of the section, we will work in the Hilbert space 

■H, :=L.2m. 

Definition of the opeiutor on the thickened, set. We want to consider the Dirichlet 
Laplaciaii on U,. Its (quadratic form is ||(iu||^. , u G 'W'(IJ,) (we could also allow 
other boundary conditions at dl x F). The transformed quadratic form is given 
by 

h,{u) := \m'u\\l = \\d{(l + eyK)-^l\)\\l^, u e ft^Ue) = fl^U,) 

A straightforward calculation already performed at various places (e. g. [ES89] 

or [DE95]) yields 

h,{u) = / / \ ^ \d^u\' + ^\dyu\' + K,\u\'']edydx (4.5) 
where the curvature induced potential is given by 

^-^^'W i(l + ey ^ 2(1 + eyK)-^ 4(1 + .-?/«:)*• ^ ' 

Note that A'^, < e < eo, is bounded from below by a constant — A^^, A^^ > 
depending only on the supremum of and £o- Using the first estimate 

in (2.7) we see that 

q^(u):=h,(u)-^\\u\\l (4.7) 
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with S = s^l'^ , Cauchy-Schwarz for \N(dxU{x, •)p < ||9a:M(a:, ■)||^ and Lemma 3.3, 
we obtain 

go(Je«)-g.(«)<0(£i/^)||u||J,i. (3.11) 
Applying Lemma 2.2 again (with A = 0), we obtain 

Aa.(Q.) = A,(e) - ^ > Afc(O) - r)A.. (3.12) 

Here, % = 0(£^/^) using (3.9), (3.11) and the upper estimate \}JyQe) < ^^^(0) + 
0(e) = 0(1) from (3.5). 

4. Curved edges 

Let us now consider a curved quantiun wave guide embedded in M- (more 
general embeddings can be treated similarly). Such spaces have already been 
analysed e.g. in [ES89] or [DE95]. We only consider a single edge here since 
one can easily replace the edge estimates without curvature by the appropriate 
estimates with curvature in the previous section^ (cf. Remark 4.1 for the precise 
assumptions on the curvature). The convergence of the discrete spectnuu of an 
infinite ciu'ved quantum wave guide has ah'eady been estabhslied in [DE95] usiug 
perturbation arguments and an asymptotic expansion (cf. also [KaP88] wliere the 
asymptotic of the first Dirichlet eigenvalue of a e-neighbourhood of a finite length 
curve in was treated). Here, in contrast, we use the variational arguments of 
Lemma 2.2 which are somehow simpler (the price being a weaker result). 

Definition of the approximating space. Suppose that 7: / — » is a smooth 
curve (e.g. is enough) with bounded derivatives parametrised by arcleugtli 
(i.e. the tangent vector 7(1-) has unit leugtli for all x 6 /). Suppose that either 
7 is a closed curve (/ = S^) or has two ends (/ = (0, 1)). 

We introduce the e-neighbourhood Ue of the curve given as the image of the 
parametrisation 

■(S-.IxF — >C/, cM^ 

{x,y) I — >j(x) + eyn{x) 

where n{x) :— (72(2;), — 7i(.x)) is orthogonal to the tangent vector 7(0;) and F = 
(— 1, 1). Define the signed curvature by 

K := 7i72 - 7271 (4.2) 

and suppose that < e < £0 := l/(2||'«||oo) where ||k||oo denotes the supremum 
of |K(a;)|, X e I. We assume in addition that 4' is a diffeomorphism. 

^More precisely: one has to show the estimates of this section for the metric 

fff = (1 - eaf(l + eyK{x)fdx^ + e'^di/ 

instead of llie melric defined in (4.3) in order 10 l alce inio acconnt. Ihe sliortcned edges due to 
the embedding. To Iteep the notation manageable we omit this fact here. 
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3. Convergence of the eigenvalues: small vertex neighbourhoods 

In this section we consider a graph with straight edges approximated by an 
open set as defined in the previous section (the case of curved edges will be 
treated in the next section). We apply the abstract comparison result Lemma 2.2 
to our concrete problem in order to show an upper and a lower bound on \k{Qe) = 

A,(A° J - Ai/e^ = A,(e) - X,/s\ 

Upper bound. We define the linear map Jq; domgo — > domg^ transmitting 
(eigen-)functions on the graph to functions on by 

W){z)-e ^^^^ (3.1) 

where ip is the first normalised Dirichlet eigenfunction on the transversal direction 
F, i.e., 

AJ^!^ = Ai(p. 

Note that flgi^ vanishes and therefore Jq/ € domq^ = '}i}{M^). We begin with 
the verification of (2.8) and (2.9). We have 



\\ml=sY."j \f(x)fdx = 0(e)\\f\\l (3.2) 

3€J '''' 



since \\<p\\f = 1. Furthermore, 
qeW) - 9o(/) 

= Y\t^ I \f\"<^^+-J I M-\imdymi-a,e)dx\ 

(3.3) 

Since If is the eigenfunction with eigenvalue Ai the latter integral vanishes and 
therefore 

QeiJof) - qoif) = E I^Till/'lll = 0{e)qo{f)- (3.4) 

Applying Lemma 2.2 with A = we obtain 

A.(e)-^<A,(0) + O(£). (3.5) 

Lower bound. For the lower bound we have to work a little bit harder. We define 
Je : dom — > dom 50 by 

{Jsu)j{x) := £™/^(Afu(a;) - p{x)Nu{x°)) (3.6) 

where 

Nu{x) ~ {u(x,-),(p) = j u{x,y)Tp(y)dy (3.7) 
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is the expectation value of u{x, •) e L2(F) corresponding to the lowest transversal 
eigenfunction ip. Here, x" depends on x and denotes the left resp. right endpoint 
of Ij if X is in the left resp. right half of Ij. Furthermore, p is a smooth function 
with < p(.t) < 1, p{x) = near the mid point of Ij and p{x) = 1 near the 
boimdary of Abusing the notation a little bit, also represents an element 
of dij. Since J^u(x") = 0, we have J^u 6 dom go- 
Again, we begin with the verification of (2.8). First, we show the following 
estimate on higher transversal modes. 

Lemma 3.1. We have 

\\vr-\{v,'P)\'<Y^{\\dvr-x4vr) 

A2 — Ai 

for V e n^(F) where A, are the Dirichlet eigenvalues of F. 

Proof. Since v — {v, if)ip is the projection onto ip-'-, the min-max principle imphes 

hr - \{v,^)\' = \\v - (v,^M' < hd{v-{v,^}^)r 

A2 A2 A2 

Since F is connected, A1/A2 < 1 and we can bring the last difference on the LHS, 
divide by (1 — A1/A2) and obtain the desired estimate. □ 

The next lemma shows that under our main assumption, eigenfunctions do not 
concentrate at the vertex neighbourhoods: 

Lemma 3.2. Assume (1.2) then 

IHlL^ Am)-A, Nl^-.--^ll"llU 
for all u e fl^M,) n n\V,^u)- 

Proof. Using the second estimate in (2.7) and the scaling of the metric (2.3) we 
have 

_2 g2 |- A 1 A 

By our main assumption (1.2), Ai/Af^(14) < 1 and the result follows as before. 

□ 

Finally, we need the following lemma. 
Lemma 3.3. We have 

e'"|Af«(x°)P<0(£)[|M«||^,^-^||«tiy 
for all u 6 ^{M^) n H^iVc^k) and 6 dIj, if ej emanates from Vt- 



BRANCHED QUANTUil WAVE GUIDES 11 

Proof. A standard Sobolev embedding theorcun gives 

INuix")]' <jHx\y)?dy < c,[\\dury^ + 

for some constant ci > (note that F = djVk). Now by the scaling of the metric 
on Vk 

\\du\\l, + \Hk=e-[e{\\du\\l^^ - ^MIJ + ]{l + X^)\\u\\'^ 
and the result follows from the proceeding lemma. □ 
Now, we want to consider the norm difference 

h\\l-\\Jsu\\l 

= Ell«ll'v.,. + E[l|-"llL- / Wu{x) - p{x)Nu{x'')\'e'-dx\. 
keK jeJ ■^'i 

The first sum can be estimated by 0{e'^)qe(u) using Lemma 3.2. For the second, 
we use 

(o + bf > (1 - 5)o' - 5>Q (3.8) 



and obtain as upper bound 

Ml.j-i'^-i) \Nu{x)fs"'dx + ^ \pix)\^dx max \Nuix°)f 

< I \\Hx,-)f -\{u{x,-),v)\^]e"^dx 
Jij 

+ - a,£) ||«||?, , + f 11,111 ma.^ \Nu{x^)f 

using Cauchy-Schwarz. Applying Lemma 3.1, the scahng of the metric on F 
in (2.3), Lemma 3.3 and setting 5 = e^l'^, we end up with the estimate 

\\ut-\\J,ut<0{e^l^)\\u\t,r. (3.9) 

For the quadratic form difference we have 

go(^£«i) - <le{v) 

< ^ e"^\\N{a,u) - p-Nuix'^l - I WdMx, ■We'^dx 

where the terms of order has been estimated with (2.7). Using 

(o + 6)^ < (1 + 5)a^ + ?^)^ 0<5<1, (3.10) 





